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ABSTRACT 

We present new optical observations of the HH 111 Herbig-Haro jet using the Gemini Multi Object 
Spectrograph in its Integral Field Unit mode. Eight fields of 5" x 3.5" have been positioned along 
and across the HH 111 jet, covering the spatial region from knot E to L in HH 111 (namely, knots E, 
E, G, H, J, K and L). We present images and velocity channel maps for the [O I] 6300+6360, Ha, [N 
H] 6548+6583 and [S H] 6716+6730 lines, as well as for the [S H]6716/6730 line ratio. We find that 
the HH 111 jet has an inner region with lower excitation and higher radial velocity, surrounded by a 
broader region of higher excitation and lower radial velocity. Also, we find higher electron densities 
at lower radial velocities. These results imply that the HH 111 jet has a fast, axial region with lower 
velocity shocks surrounded by a lower velocity sheath with higher velocity shocks. 

Subject headings: ISM: Herbig-Haro objects — ISM: jets and outflows — ISM: kinematics and dy¬ 
namics — ISM: individual (HH 111) 


1. INTRODUCTION 

HH 111 (discovered by Reipurth 1989) is one of the two 
most remarkable and better studied Herbig-Haro (HH) 
jets, the other one being HH 34 (shown to be a jet by 
Reipurth et al. 1986). The past observational studies of 
HH 111 include: 

• ground based (Reipurth et al. 1992; Podio et al. 
2006) and HST (Raga et al. 2002a) high and low 
resolution long-slit spectra, 

• optical (ground based: Reipurth et al. 1992; HST: 
Hartigan et al. 2001) and IR (Goppin et al. 1998) 
proper motions, 

• evidence of multiplicity of the outflow source (Gre- 
del & Reipurth 1993; Reipurth et al. 1999; 
Reipurth et al. 2000; Noriega-Crespo et al. 2011), 

• discovery of an associated giant jet (Reipurth et al. 
1997a), 

• detection of an associated, very well collimated 
molecular outflow (Cernicharo & Reipurth 1996; 
Nagar et al. 1997; Hatchell et al. 1999; Lefloch 
et al. 2007) 

Numerical simulations of variable jets calculated 
specifically for modelling the observational properties of 
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HH 111 were presented by Masciadri et al. (2002) and 
Raga et al. (2002b). 

Optically, the HH 111 system has a one-sided jet that 
appears at ^ 15" from an obscured source (detected at 
radio wavelengths by Reipurth et al. 1999), extending 
to a distance of ^ 40" W from the source. However, a 
much more symmetric jet/counterjet structure extending 
down to the position of the source is observed in ground 
based (Gredel & Reipurth 1994; Davis et al. 1994), HST 
(Reipurth et al. 1999) and Spitzer (Noriega-Crespo et al. 
2011) IR images. At angular scales > 2' from the source, 
both of the outflow lobes are detected optically, with 
total extent of ^ 2° for the whole system (see Reipurth 
et al. 1997a). 

The kinematics of the HH 34 jet have been studied 
spectroscopically (for a limited set of optical emission 
lines) with full spatial coverage by Beck et al. (2007) and 
Rodriguez-Gonzalez et al. (2012). Our present paper 
describes similar observations (i.e., high resolution spec¬ 
troscopy with full spatial coverage), but for the HH 111 
jet. In our Integral Eield Unit (lEU) Gemini South spec¬ 
tra the [S H] 6730, 6716; Ha; [N H] 6583, 6548 and [O I] 
6360, 6300 are detected, and have a high enough signal- 
to-noise so that velocity channel maps can be generated. 

The paper is organized as follows. In section 2 we 
describe the observations and the data reduction. The 
results are described in section 3. Einally, the work is 
summarized in section 4. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. The observations 

HH 111 was observed at GEMINI North observatory 
on Oct 2007 and Nov 2010, using the GMOS instrument 
in lEU mode (Allington-Smith et al. 2002) under the 
program GN-2007B-Q-9. We used the lEU in the single 
slit mode with the R831_G5302 grating, giving a R = 
4396 resolution at 7570A. In the single slit mode, the 
lEU field of view is 5" x 3"5, each lens covering 0"2 on 
the sky. 
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TABLE 1 

Target coordinates and offsets 


Target 

RA 

(h:m:s) 

DEC 

Ap 

(") 

Aq^ 

(") 

IRAS 05491+0247‘> 

05:51:46.25 

02:48:29.5 

- 

- 

field 1 

05:51:43.7 

02:48:34.1 

- 

- 

field 2 



0 

4.4 

field 3 



0 

8.8 

field 4 



2.9 

5.6 

field 5 



-2.9 

5.6 

field 6 



0 

-4.4 

field 7 



0.214 

-9.07 

field 8 



0.214 

-13.43 


^The offsets Ap and Aq are given in arcsec, and refers to dis¬ 
placements perpendicular (p) and parallel (q) to the HH 111 PA; 
here defined as ~ 83°. Field 1 is the reference to the displacement 
shifts. 

^Coordinates of the HH 111 IRS (Reipurth et al. 2000). 



Fig. 1.— Sketch of the eight fields observed with the GMOS- 
IFU superimposed on the Ha pre-image of HH 111 taken with 
the Gemini North Telescope. The observed fields are labeled by 
numbers from 1 (centered on the brilliant knot J) to 8 (finishing the 
mosaic in the double knot E). The knots inside the HH 111 jet are 
labeled by capital letters from E to L, following the nomenclature 
provided by Reipurth (1989) and Raga et al. (2002a). The offsets 
between each field are indicated in Table ^ The N-E axes are 
indicated in the Figure, as well as a distance scale. 


The observations has made under exceptional seeing 
conditions. Using the R-band images, the seeing esti¬ 
mates range from 0"44 to 0"51 for the first epoch obser¬ 
vations (2007) and from 0"55 to 0"60 for the second one 
( 2010 ). 

Ha narrow-band filter pre-images (on and off) were 
also taken (on August 22nd, 2007) to position the IFU 
field with higher accuracy. In order to cover the tar¬ 
get, we observed the 8 fields shown in Figure Each 



Fig. 2.— Left: Ha sub-cube of field 8 without any filtering process 
(3"5 X 5"0). Middle: The sum of the first two eigenvectors, after 
PGA analysis of the data cube in the spectral region defined by Ha, 
but excluding the emission. Right: First image (left panel) minus 
the low-frequency noise (middle panel). The results show that the 
instrumental fingerprint has been removed from the original data, 
without any introduction of artificial structures. 

field was observed with a total exposure of 1200s (three 
400s exposures, which were then median averaged to re¬ 
duce the cosmic ray contamination). Fields 1 to 6 were 
observed on October 19th 2007 and fields 7 and 8 on 
November 16th 2007. Fields 1 to 6 were re-observed on 
November 2nd 2010, due to bad CCD settings during 
the original observations. Overlaps between the different 
fields have been set to 0"4 (corresponding to two rows of 
IFU lenses). 

Table gives the coordinates of the HH 111 IRS 
outflow source (IRAS 05491+0247; see, Rodriguez & 
Reipurth 1994; Rodriguez et al. 1998; Reipurth et al. 
2000) and the exact offsets of the GMOS-IFU observed 
fields. The offsets are in arcsec and are taken with re¬ 
ject to the center of field 1 (center of knot J, see Figure 
Q. The offsets are perpendicular (Ap) and parallel (Ag) 
to the HH 111 outflow axis, which is at a PA=—83.94° 
position angle (Raga et al. 2002a). 

2.2. Data reduction 

Data have been reduced using special reduction pack¬ 
age provided by Gemini Staff using the IRAF reduc¬ 
tion package R All raw images have been bias corrected, 
trimmed and nat fielded. Flatfield images have also been 
used to located the positions of the 1500 lenses on the 
frame. Twilight images were used to estimate the grat¬ 
ing response and after the extraction of spectra for each 
lens, each one has been corrected using this normalized 
response. Arcs, from the GuAr lamp, have also been 
taken for the wavelength calibration. Using the bright 
01 night sky line at 5577.338 A, we have estimated the 
wavelength calibration accuracy to 0.1 A. The last step 
of the reduction was the sky subtraction using the field 
located at V from the science field. Finally the data 
cube has been created with the GFGUBE task using a 
spatial resampling of 0.1" per pixel. The total spectral 
coverage goes from 4835.896 A to 6957.800 A. The spec¬ 
tral sampling is 0.339 A per spectral pixel. We then 
extracted sub-cubes near the Ha, [O I]AA 6300,6364, [N 
H]AA 6548,6584 and [S H]AA 6716,6731 emission lines, for 
each the eight observed fields. 

Eurther image treatment has been conducted using the 

^ http://www.gemini.edu/node/10795 

^ Image Reduction and Analysis Facility is a software developed 
by the National Optical Astronomy Observatory - iraf.noao.edu 
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Fig. 3.— From left to right: Mosaics of [S II]A6716+A6731 for the integrated line profile (left), and for three different radial velocities 
intervals: high velocity (HVC, medium velocity (MVC) and low velocity (LVC) components. The intervals are defined by: Vrad < —100 
km s“^ for HVC, —100 < v^ad < “"^0 (in km s“^) for the MVC and v^ad > “"^0 km s“^ for the LVC. The coordinate axes are in arcsecs 
(the origin coincides with the position of the driving source) and the colorbar are in arbitrary units. 



d (arcsec) d farcsec) 

Fig. 4.— The same as in Figure for Ha. 

pipeline developed by Menezes et al. (2014; see also 
Ricci, Steiner & Menezes 2011). It consist in i) a spatial 
filtering process to remove high-spatial frequency noise, 
ii) a Principal Component Analysis (PCA) to remove 
instrumental spatial fingerprints and iii) a Richardson- 
Lucy deconvolution process. In order to remove the high 
spatial frequency noise, we first calculate discrete Fourier 
transforms of the images resulting from each data cube, 
and apply a low-pass band filter in the frequency space 




d [arcsecl d farcsec) 


(Gonzalez & Woods 2002). To this effect we follow the 
procedure described in Menezes, Steiner & Ricci (2014). 
In particular, we use a Butterworth H{u^v) filter of or¬ 
der n = 6 (see equation 7 in Menezes et al. 2014). 
An inverse Fourier transform is then applied to the fil¬ 
tered data. We should mention that similar noise filter¬ 
ing process has been already applied in the context of 
long slit spectroscopy of HH jets (Raga & Mateo 1998). 
The high-frequency cleaned data still show low-frequency 
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Fig. 5.— The same as in Figure]^ but for [N IIJA6548 + A6583. 




d (arcsec) d fancsec) 


Fig. 6.— The same as in Figure[^but for [O IJA6300 + A6360. 

noise which we remove using the Principal Component 
Analysis (PCA; see Ricci, Steiner & Menezes 2011 and 
Steiner et al. 2009 for a detailed discussion and method 
presentation). The scope of this paper is not to present 
this technique. However, it is worth to mention briefly 
how it works. The PC A technique allows us to describe 
the data cube as a linear combination of an orthogonal 
basis, which are the eigenvectors of a 2D covariance ma¬ 
trix. Such a covariance matrix is obtained after scaling 
all the intensities in the 3D data cube with the transfor¬ 




d (arcsec) d farcsec) 




d (arcsec) d fancsec) 


mation: 


/3 = - 1) + i (1) 

in a data cube with the spatial dimensions n = /i x z/: 
l<i</i, l<j<z^. For each spaxel intensity, 
lijx^ one defines a 2D intensity matrix which will 
be the subject of the PCA analysis. Its covariance ma¬ 
trix is calculated, and the eigenvectors with their asso¬ 
ciated variances are obtained (e.g., Steiner et al. 2009). 
The higher its variance, the more representative of the 
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Fig. 7. — Low-velocity and high-velocity components mosaic sub¬ 
traction in Ho: emission. Left: LVC - 2xHVC. Riqht: LVC - 
4xHVC. 

data is a given eigenvector. As the orthogonality is en¬ 
sured, we can compose a tomogram, which is an image 
that represents the projection of the original data onto 
a selected eigenvector (or eigenvectors). In Figure!^ we 
show a specific observed field (field 8 in Figure fT]h for 
a given emission line (Ha), before any filtering process 
(a), after the filtering process but extracting the emis¬ 
sion line from data (b), and the (a)-(b) subtraction in 
(c). We should mention that the tomogram in (b) was 
obtained by adding the first two eigenvectors (without 
the emission line), which we assume to correspond to the 
low-frequency noise, or the instrumental fingerprint, as 
described in Ricci et al.(2011). 

Finally, we have applied a Richardson-Lucy decon¬ 
volution procedure with a Gaussian PSF, with a con¬ 
stant FWHM of 0"44 for fields 7 and 8, and a constante 
FWHM of 0"55 for the fields 1 to 6, as suggested by 
the seeing in our observations. Six interactions has been 
used in order to deconvolve the original data with the 
PSF (see equation 10 in Menezes et al. 2014). The 
whole procedure (spatial re-sampling, low-noise filtering 
process, instrumental fingerprints removal technique and 
Richardson-Lucy deconvolution) is explained in detail in 
Menezes et al. (2014), who has developed the techniqu^ 

3. RESULTS 

3.1. Images and velocity channel maps 

We present the spatially resolved line profiles as “ve¬ 
locity channel maps” in three velocity intervals: 

• HVC: a high (negative) velocity channel, with ra¬ 
dial velocitie^ < —100 km s“^, 

• MVC: a medium velocity channel, with —100 km 
S~^ < Vrad < “70 km s“^, 

• LVC: a low velocity channel, with v^ad > ~70 km 
s“^. 

® The routines are available electronically at the following ad¬ 
dress: http://WWW.astro.iag.usp.br/$\sim$PCAtomography/ 

^ The velocities has not been corrected for the LSR velocity of 
+23 km s“^ Reipurth 1989. 


We also present an image (for each spectral line) which 
consists of an addition of these three velocity channel 
maps. The channel maps and the images have been com¬ 
puted for mosaics composed of the 8 observed IFU fields 
(see Figure 1), and cover the emission of the HH 111 jet 
at distances ~ 24" ^ 50" from the outflow source. The 
mosaics have been built taking into account the offsets 
(given in Tablel^ 

In Figures [3[^ and |^we show the images and 
velocity channel maps for ^ HJA6716+A6731, Ha, [N 
HJA6548+A6584 and [O I]A6300-fA6360^espectively. In 
the [S H] image (left frame of figure 1^ we see knots 
E through L, which can be clearly recognized compar¬ 
ing the image with the identifications of Reipurth et al. 
(1992) and Reipurth et al. (1997b). 

In the channel maps of all of the observed lines, it is 
clear that there is a general trend of increasing jet width 
as a function of decreasing velocity (i.e., in all of the ob¬ 
served emission lines the jet becomes narrower as we go 
from the HVC to the MVC and to the LVC, see Figures 
3][6 ). This result is in qualitative agreement with the ob¬ 


servations of Riera et al. (2001), who found velocities 
decreasing away from the outflow axis in two long-slit 
spectra cutting the HH 111 jet at the approximate po¬ 
sitions of knots D (not detected in the present observa¬ 
tions) and F. 


3.2. Width vs. radial velocity dependence 

In order to illustrate the broadening of the jet at lower 
(i.e., less negative) radial velocities, in Figure 7 we show 
the LVC—HVC channel map subtraction for the Ha line. 
This subtraction map shows that at all positions along 
the observed region of the HH 111, the high radial veloc¬ 
ity emission is restricted to a narrow, central region of 
the jet cross section. 

One can quantify the width vs. radial velocity depen¬ 
dence by measuring the FWHM of the emission on the 
HVC, MVC and LVC channel maps. For all of the emis¬ 
sion lines, we measure the widths on the three channel 
maps at the positions of knots E, E, G, H, J and L, in¬ 
tegrating the emission in bins of 1" along the jet bin 
(centered at the positions of the knots). The results of 
this procedure are shown in Eigure 8, which shows the 
following feature^ 

• there is a systematic increase of width as a function 
of decreasing radial velocity for all the knots, and 
in all emission lines, 

• with the exception of the very broad knot E, there 
is a general trend of increasing width vs. distance 
from the source, 

• in the lower excitation lines ([O I] and [S H]), the 
widths are somewhat lower than in the higher ex¬ 
citation lines (Ha and [N H]). 


This last effect has been noted by Reipurth et al. 
(1997), who note that the Ha emission is broader than 

® The determination of the FWHM for the knots Gl, J, K1 and 
L, using the [N II] and [O L lines show high uncertainties due to 
their low S/N (see I n ures [m and. In a few cases (knots GI, H 
and KI; see Figure]^ we were not able to determine the [N II] 
FWHM for some (if not all) velocity components. 
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Fig. 8. — FWHM (in arcsecs) of the emission on the HVC, MVC and LVC channel maps for all of the emission lines (indicated on the 
top of the figure), at the positions of knots E, F, G, H, J and L (indicated on the right of the figure). The FWHM has been obtained 
integrating the emission in bins of along the jet (centered at the positions of the knots) and fitting a Gaussian to each knot longitudinal 
profile. For [N H] line we found it difficult to adjust a profile due to low S/N in some knots and for some velocity components (see also 
Figure!^. For these cases, the FWHM appears in the histogram as an zero. The Ho: and [N H] histograms for the knot L has the FWHM 
axis divided by a factor of 2. 


the [S II] emission (see their Figure 6, in which they 
present an H(a-[S II] subtraction map). In our obser¬ 
vations, we see that in general the [N II] emission has 
widths comparable to the Ra widths, and that the [O I] 
emission has narrower widths, comparable to the ones of 
the [S II] emission. 

3.3. The [S II] 6716/6731 line ratios 

In Figure we present the spatial distribution of the [S 
II]A6716/[S n]A6731 line ratio for the integrated line pro¬ 
files and for the HVC, MVC and LVC maps. This ratio 
gives the electron density {ue) of the medium Osterbrock 
(1989). For calculating rzg we assume a temperature of 
10^ K, consistent with the temperature of ^ 1.3 x 10^ K 
found by Podio et al. (2006) for the HH 111 jet. 

Superimposed on the spatial distribution of the elec¬ 
tron density are the (logarithmic spaced) iso-contours of 
the [S II]A6716 + A6731 integrated intensity. The inte¬ 
grated map (left) shows that the electron density ranges 
from ^ 1.5 X 10^ cm“^ at knot E (at ^ 26 arcsec from the 


outflow source) up to 5 x 10^ cm in knot F (at ^ 2"5 
from knot E; see Eigurej^. It then drops to ^ 2 x 10^ 
cm“^ between knots E and H. In knot H, the electron 
density has a strong peak of ne ^ 3.5 x 10^ cm“^. 

Interestingly, the electron density maps of the three ve¬ 
locity channels show systematically increasing densities 
for decreasing radial velocities. This is can be clearly 
seen at the positions of knots E, E and J. Also, an in¬ 
crease in Uq can be seen in the region between knots E 
and G when going from the MVC to the LVC electron 
density maps (see the two frames on the right of Eigure 
9). 

3.4. Knot E 

We now focus on the interesting structure of knot E (at 

26" from the source, see Eigures 3-7). The structures 
observed in the images and channel maps of this knot are 
shown in Eigure 10. 

We see that in the Ra and [N II] images (2nd and 3rd 
rows of the far left column of Eigure 10) knot E shows 
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Fig. 9. — From left to right: Mosaics of [S II]A6716/[S IIJA6731 line ratio for the integrated line profile, and for three different radial 
velocities intervals: high velocity (HVC), medium velocity (MVC) and low velocity (LVC) components. The intervals are defined by: 
'^rad < —100 km s“^ for HVC, —100 < v^ad < “70 (in km s“^) for the MVC and Vrad > “70 km s“^ for the LVC. The coordinate axes 
are in arcsecs and the colorbar indicates the electron density in cm“^. 


two sidc-by-sidc peaks. These peaks have been observed 
in previous ground based (Reipurth et ah 1992) and HST 
(Reipurth et ah 1997) Ra images of HH 111. Most inter¬ 
estingly, the two side-by-side peaks are well separated in 
the LVC Ra and [N H] maps (right column of Figure 10), 
approach each other in the MVC map, and merge into a 
single, central peak in the HVC map. Therefore, we find 
that the two side-by-side peaks (previously observed in 
Ra images, see Reipurth et al. 1997) actually enclose a 
fainter, higher radial velocity emission region. 

In the lower excitation [S H] and [O I] lines, knot E 
shows two low intensity side-by-side peaks only in the 
LVC maps (1st and 4th rows of the far right column of 
Figure 10), and these peaks are not present in the MVC 
and HVC maps. The images of knot E in these lines (1st 
and 4th rows of the far left column of Eigure 10) show a 
conical structure which resembles the MVC emission. 

4. SUMMARY 

We have obtained lEU spectra of the HH 111 jet at 
distances of 24 ^ 49" from the outflow source. Erom the 
spectra, we obtain position-velocity cubes for the [S H] 
6716, 6731; Ra; [N H] 6548, 6583 and [O I] 6300, 6360 
emission lines. 

We find a number of interesting results: 

• we find that for all emission lines we have an in¬ 
crease in the observed width of the jet for decreas¬ 
ing radial velocities. This effect can be seen directly 
in the velocity channel maps (Eigures 3-6), in sub¬ 
tractions of pairs of velocity channel maps (Eigure 
7) or more quantitatively in the widths measured 
for the jet knots (Eigure 8). A broadening in the 
HH 111 jet for decreasing velocities has been previ¬ 
ously observed by Riera et al. (2001), and a qual¬ 
itatively similar effect is seen in the CO emission 
(but at widths in excess of ^ 10", see Lefloch et al. 
2007), 


• we find that the HH 111 jet is narrower in the 
lower ([S H] and [O I]) than in the higher {Ra and 
[N H]) excitation lines at the vicinity of some knots 
(namely, knots E, E K1 and L). This effect has been 
previously seen in the Ra and [S H] HST images 
of Reipurth et al. (1997). These previous observa¬ 
tions suggested that we might be seeing a sheath of 
non-radiative shocks (seen only in Ha) surrounding 
the HH 111 jet beam, but the detection of broad [N 
H] emission indicates that the shocks in the sheath 
are indeed radiative, 

• from the [S H] 6716/6731 line ratio we systemat¬ 
ically find larger electron densities at lower radial 
velocities (see Eigure 9), 

• the “twin peak” structure of knot E is clearly seen 
in the Ha and [N H] images, but these two peaks 
merge into a central peak in the [O I] and [S H] 
images. In all of the lines, the twin peaks are seen 
at low radial velocities (see the LVC maps of Eigure 
10), but they approach and merge into a single peak 
at higher radial velocities. 

These characteristics can be qualitatively explained in 
terms of an “internal working surface” model, in which 
the knots are produced as a result of an outflow velocity 
time-variability. In such a model, the variability pro¬ 
duces pairs of shocks that travel down the jet beam, pro¬ 
ducing a high pressure region that ejects material into the 
cocoon of the jet. The observations of larger widths for 
the higher excitation emission would then be interpreted 
as higher velocity shocks produced by the material in this 
sideways ejection against a slower moving cocoon. The 
low excitation emission would be interpreted as coming 
from lower velocity shocks within the jet beam (produced 
by the outflow velocity variability). 

The shock structure of internal working surfaces was 
described by Raga et al. (1990), and numerical Simula- 
































Cerqueira et al. 


-10 12 
d (arcsec) 


-10 12 
d (arcsec) 











Fig. 10.— From left to right: Integrated emission, HVC, MVC and LVC for (from top to bottom) [S II] , Ha, [N II] and [O I] for field 8 
which contains the knot E. The axes are in arcsec and indicates the distance from the driving source. The maps are in arbitrary units (as 
seen at the colorbars). 
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tions have been done for reproducing Ha images (Raga et 
al. 2002b) and long-slit spectra (Masciadri et al. 2002) of 
the HH 111 jet. The data presented here clearly provide 
a substantial increase in the observational constraints, 
which should be a challenging test for future models of 
jet knots produced by an outflow variability or by other 
mechanisms (see e.g., Micono et al. 1998). In particular, 
our data will give interesting constraints on the nature 
of the cross section of the outflow (Raga et al. 2011) 
or of the possible presence of two distinct components 
(Tesileanu et al. 2014). 

We can speculate that MHD effects might be present, 
and might help to shape the morphology of the HH 111 
emission knots as well as their line profiles. It is worth 
to note that recent and sophisticated MHD models pre¬ 
dict that beam recollimation (by the hoop stress) plays 
an important role in order to enhance the brightness 
near the jet axis, as well as to accelerate the inner jet 
beam, making it faster than the surrounding jet mate¬ 
rial. This general result was obtained by Tesileanu et al 
(2014), Hansen, Frank & Hartigan (2014) and Staff et 
al. (2015), in spite of the differences in the initial setup 
in these works. In Hansen, Frank & Hartigan (2014), a 
pure toroidal magnetic field is used to study the jet prop¬ 
agation, while the disk wind solution is incorporated in 
the models from Staff et al (2014) and Tesileanu et al. 
(2014), who actually introduced a two component (stellar 
jet plus disk wind) flow. 

Another feature shown by some HH jets (including 
some of the knots studied here) is an increase of the 
FWHM of the jet beam with distance from the driving 
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source. This feature has been predicted in the models 
presented by Tesileanu et al. (2014) and compared with 
available data from RW Aur, HH 30 and HL Tan, while 
the solutions presented in Staff et al. (2014) were com¬ 
pared with data from DG Tan, HN Tan, RW Aur and 
UZ Tan. In both cases, the distances covered by the 
simulations were suitable for the study of the microjets 
(^ 100 AU). Whether or not these effects are still impor¬ 
tant farther away from the driving source (the knot E is 
at ^ 10 000 AU from the VLA 1 source) is an issue that 
should be addressed properly in future work. 
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